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Abstract 



We study the off-equilibrium dynamics of a particle in a general N- 
dimensional random potential when N — ► oo. We demonstrate the existence of 
two asymptotic time regimes: i. stationary dynamics, ii. slow aging dynamics 
with violation of equilibrium theorems. We derive the equations obeyed by 
the slowly varying part of the two-times correlation and response functions 
and obtain an analytical solution of these equations. For short-range corre- 
lated potentials we find that: i. the scaling function is non analytic at similar 
times and this behaviour crosses over to ultrametricity when the correlations 
become long range, ii. aging dynamics persists in the limit of zero confining 
mass with universal features for widely separated times. We compare with the 
numerical solution to the dynamical equations and generalize the dynamical 
equations to finite N by extending the variational method to the dynamics. 
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I. INTRODUCTION 



A major theoretical challenge, as well as an important issue for numerous experimental 
systems, is to understand the non-equilibrium dynamics of elastic manifolds, with or without 
internal periodic structure, in quenched random media 1 . The corresponding models, such as 
the model of an interface subjected to quenched impurities or the Sine Gordon model with 
phase randomness, have been studied for some time. Some progress has been made in the 
description of the statics 2-6 but little is still known about the dynamics, especially about the 
non-equilibrium features of the relaxations. It is by now well-established that these systems 
exhibit glassy behaviour which presents similarities with the glassiness of spin-glasses, such 
as slow relaxations, history dependence and strong sensitivity to changes in the external 
parameters 7 ' 8 . 

Recently progress has been made in solving the non-equilibrium dynamics of mean- 
field models of spin-glasses 9-13 . A new method to study analytically the large-times off- 
equilibrium dynamics has been introduced 10 ' 11 . This purely dynamical method has on the 
one hand some formal connections with the replica symmetry breaking schemes introduced 
to study the static properties, i.e. the Gibbs measure, of these models. On the other hand, 
the method yields interesting additional information specific to the dynamics. It allows 
us to find a solution which exhibits two asymptotic time regimes: a stationary dynamics 
regime for large but similar times, and a slow aging regime for widely separated times. 
Most importantly it then allows to establish contact with the essentially out of equilibrium 
experimental observations, namely the slow relaxations and the aging effects 7 . It is then 
natural to extend this method to study the dynamics of a broader class of systems with slow 
relaxations including the random manifold problem. 

For problems such as the manifold in a random medium, the mean-field limit corresponds 
to the dimension of the embedding space N going to infinity, i.e. it is represented by a field- 
theory with a large number of components. In this case one can derive a closed set of 
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dynamical equations and the above dynamical method can be applied yielding exact results. 
Realistic systems, however, are embedded in a finite dimensional space, i.e. N is finite, 
and one must make some approximation to obtain closed dynamical equations. A way to 
obtain these equations is to extend the Gaussian Variational Approximation (GVA) to the 
dynamics. 

The model of a manifold of internal dimension D embedded in a random medium 
of dimension N is described, in terms of a N component displacement field </>,</> = 

(<f>l,<f>2,---,<f>N) 



by the Hamiltonian 4 ' 5 



H= I d D x 



XV 



fi is a mass, which effectively constraints the manifold to fluctuate in a restricted volume of 
the embedding space. V is a Gaussian random potential with correlations 

V(<l>, x)V (</>', x') = -N6 D {x - x') V ^ • (1-2) 

Our aim is to study the out of equilibrium dynamics of this model for a general random 
potential. The dynamics that we consider is the Langevin dynamics 

dd>(x,t) 8H . . 

^f J = -^7) + "(- i » < L3 > 

with (i,„(:c, t),i s (x', (')) = 2T 6„g 6 D (x - x')6(t - f). 

To summarize, in this paper we concentrate on the problem of a particle moving in a N 
dimensional random potential. It corresponds to the limit of a manifold with zero internal 
dimension 14 ' 15 D = 0. The zero dimensional case is a necessary step before studying the 
model in finite D and a good trial ground for the method. Schematically we do the following. 
Firstly we study analytically the exact dynamical equations in the mean-field limit N — > oo. 
Secondly, we derive the dynamical equations for finite N within the extended GVA and we 
apply the same dynamical method to the approximated set of equations. We shall present 
the analysis of the finite dimensional models (D > 0) with the necessary extensions of the 
method in a separate publication 16 . 



The quantities of interest in the long-time off-equilibrium dynamics are the two time 
correlation and response functions. In the D = model we define 



^ V" ) f=o 

where fit 1 ) is a small perturbation applied at time t'. For the models we consider here 
it is also convenient to use the mean squared displacement correlation function B(t,t') = 



Infinite dimensional models N — > oo 

There are several models which can be studied for N — > oo and D = corresponding to 
different choices of the random potential correlation V(z). A common choice is 



In what follows we shall refer to this model as the power law model. There are two physically 
distinct cases. If 7 < I the correlations grow with the distance z and the potential is called 
'long-range'. If 7 > I the correlations decay with the distance and the potential is called 
'short-range'. The statics of these models has been studied with the replica trick 5 ' 17 . The 
short range case is considered to be solved with a one step replica symmetric ansatz while 
the long-range case needs a full replica symmetry breaking scheme. 

One can also define spherically constrained models letting the mass // become a time- 
dependent function fi(t) related to the Lagrange multiplier enforcing the iV-spherical con- 
straint. Model (1 .1) then represents a manifold of internal dimension D embedded on the 
surface of a iV-dimensional sphere and it can be related to spherically constrained spin-glass 
models. 

In particular, the p-spin spherical spin-glass model defined by the Hamiltonian 18 






c(m) + c(f,f) - 2C(mo = am - m?)- 





N 



V (<f>)= J ii,-,i P <f>ii ••• < f>ip 




ii<...<i p 
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with the additional spherical constraint YliLi $ = N and taken from a Gaussian dis- 

tribution is contained in this family of models. It is recovered for D = and the correlation 
of the random potential V given by: 

V(*) = - 5 (l-2) • (1-7) 

The replica analysis yields an exact replica symmetric solution 19 when p = 2. However, this 
is a marginal case and an infinitesimal departure from p = 2 makes the exact solution be 
one step replica symmetry breaking 18 . The out of equilibrium dynamics has been studied in 
Ref. [33,13] for p = 2 and in Ref. [10] when p > 3. 



We study these models as follows 10-13 . The thermodynamic limit N — > oo is taken 
initially and then, only afterwards, the limit of large times is eventually taken. The dynamics 
starts at a finite initial time t = 0. The initial condition is chosen at random from a Gaussian 
distribution. The models we study have many ergodic components and the system is not 
able to jump from one component to another since, by definition, the barriers separating 
different ergodic components are divergent and hence unsurmountable. However, almost all 
the initial conditions do not lead the system to equilibrium at large time. A clear realisation 
of this is observed studying the dynamics of the p = 2 spherical spin-glass 13 starting from a 
generic initial condition for which the system evolves indefinitely without ever reaching any 
of the two equilibrium states. 

The two assumptions that describe quantitatively this scenario are the weak ergocity 
breaking hypothesis 20 ' 10 and the weak long-term memory hypothesis 10 ' 11 . The former states 
that even in the limit of large time the particle is always able to escape from its previous 
positions and never gets stuck in a local equilibrium. The latter states that at a sufficiently 
large time t the particle forgets any small perturbation applied during any previous finite 
time interval. For instance the 'remanent magnetisation' associated to a field applied during 
a finite time interval [0, t w ] decays to zero for a sufficiently large t. For mean-field spin-glass 
models the relaxations slow down as time passes though the systems do not reach any kind 

6 



of equilibrium. There is no equilibration time t eq such that for all the subsequent times 
the 'equilibrium-dynamics' theorems, i.e. time-translation invariance (TTI) for all the cor- 
relation functions and the fluctuation-dissipation theorems (FDT) hold. The systems do 
not visit the equilibrium states at any time. Importantly enough, this scenario is what is 
observed during experimental times in real spin-glasses. Despite being out of equilibrium, 
mean-field spin-glass models as well as real spin glasses reach an asymptotic (large times) 
regime for which some general features can be demonstrated. Generalisations of the equi- 
librium theorems to the large-times out of equilibrium relaxations have been proposed in 
Refs. [10,11] and we shall apply them here to study model (1.1) for a general V. 

It is important to note that this scenario is completely different from the one proposed by 
Sompolinsky and Zippelius 21 ' 22 to study the equilibrium properties of mean-field spin-glass 
models using a dynamical approach. In the so-called Sompolinsky-dynamics one takes a large 
but finite system (finite N). One then chooses the initial time to be minus infinity, in such 
a way that one assumes that at finite times the system has arrived at a certain equilibrium 
state. One namely takes the inverse order of limits as above: limjv^oo lim^oo to allow 
for jumping from one ergodic component to another 22 . In this approach time translational 
invariance is always assumed and hence aging effects are not captured. Moreover, the method 
has some problems since when N is finite and times diverge with respect to N quantities 
such as the the auto-correlation function become non self-averging 23 . The model (1.1)-(1.5) 
has been studied previously in this way by Kinzelbach and Horner 24 ' 25 . 

In this paper we find that, in the low temperature phase, for large times t > t' such 
that it — t')/t' << 1 the particle undergoes stationary dynamics B(t,t') = Bpit — t') where 
Bpit — t') grows from to the limiting value q. For larger time separation it — t')/t' = 0(1) 
there is an aging regime where B(t,t') further grows slowly from q to b . The behaviour in 
the aging regime depends on the range of the correlations of the potential. For long-range 
correlations we find that the one time quantities converge to the values predicted by the 
statics. This confirms the numerical observations of Ref.[12] for model (1.1) with V given 
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by (1.5) and 7 < 1. We demonstrate that B(t,t') satisfies ultrametricity in time, and we 
explicitly compute the rate of violation of the FDT theorem at large times. For short- 
range correlated potentials we find that one time quantities which involve the aging regime, 
such as the asymptotic energy density or b 0} differ from their equivalents in the statics. 
The dynamical phase diagram is different from the equilibrium one and, in particular, the 
transition in the dynamics survives in the limit of vanishing mass //—*■(). We find that in 
the aging regime B(t,t') = B[h(t')/h(t)] where generically the scaling function B[X] is non 
analytic when A — > 1 (t 1 /t — > 1) with a non trival exponent a. Furthermore in the limit 
of vanishing mass — > we show that, for all short-range models, B[X] ~ q + ln(l/A) for 
A << 1. This remarkable form implies that the response function at large time separation 
becomes a function of t' only. This is consistent with some numerical results that we present. 

The parameter 7 in the definition of V, Eq. (1.5), tunes from long-range correlated 
potentials (7 < 1) to short-range correlated potentials (7 > 1). Interestingly enough, this 
allows us to explicitly show how the ultrametric dynamical solution is approached when 
7 — > 1 from above. 

The limit of zero mass raises some fundamental questions. For fi > 0, C(t,t) reaches a 
finite limit when t — > 00. By contrast, at strictly // = 0, C(t,t) grows unboundedly with 
time. Although the corresponding unbounded diffusion process deserves further study 16 , we 
show that the method still applies and the limit //—*■() can safely be taken for quantities 
such as B(t,t') and R(t,t'). 

A first analysis of the out of equilibrium relaxation of model (1.1) at D = with long- 
range correlations given by Eq.(1.5) has been carried out by Franz and Mezard 12 . They 
have solved numerically the mean-field dynamical equations and have proposed a scenario 
of non-overlapping time domains to account for their numerical results. 
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Finite dimensional models N < oo 



The model (1.1) with finite N (and D = 1, 2, 3) has several physically interesting realiza- 
tions, such as the problem of a manifold pinned by impurities for which the correlator is of 
the form (1.5). It arises in the study of interfaces in a random field as well as glassy phases 
of vortices in high T c superconductors. Similarly the Sine Gordon model with random phase 
disorder (RSGM) arises in the study of several problems including quenched disorder. For 
example it describes the glass transition of a surface of a crystal deposited on a disordered 
two-dimensional substrate. It is also related to the vortex-free XY model in a random field. 
It is defined in term of a N = 1 phase field ip by: 



H = d 



C jJj 

- (Vtjj(x)) 2 - Ci(x) cos(V'(x)) - ( 2 (x) sm(ip(x)) + - 4> 2 



X8) 



where Ci( x )Cj( x ') = ^ij^{ x ~ x ')- This model corresponds to the choice of correlator: 



V(V>, x)V{ip', x 1 ) = -6 D (x - x 1 ) cos(V> - • (1-9) 

The statics of the RSGM has been studied using renormalization group techniques 2 , a Gaus- 
sian approximation (variational method) complemented by the replica trick 6 and with ex- 
tensive numerical simulations 3 . 

The statics of model (1.1) for general D and N has been studied using a Hartree, i.e. a 
Gaussian variational approximation 5 . The corresponding mean-field equations are identical 
to those for N infinite, apart from a replacement of the random potential correlator V(x) by 
V(x) defined as V((</>)q) = (V(</> 2 ))o where (-) denotes a Gaussian average over <f>. In this 
paper we extend the GVA to the dynamics by performing a Gaussian decoupling approxima- 
tion in the exact dynamical equations of motion. As it is detailed in Appendix A, this shows 
explicitly that the same replacement of V(x) by V(x) holds in the dynamics. For instance 
the RSGM dynamics is described by Eq. (1.5) with V(x) — > V(x) = — A exp(— x/2). 

The zero dimensional version D = has been studied as a toy model of the above 
problems 14 ' 15 . Even if one does not expect strictly speaking a glass transition in D = 0, 
this model exhibits at low temperatures several features of a glass which are present in its 
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higher dimensional versions. In the statics perturbation theory breaks down because ol the 
large number of metastable states and the importance of rare fluctuations. In the dynamics, 
there is a finite ergodic time t eq beyond which equilibrium dynamics is established for all 
subsequent times. This time however can be very large since barriers grow as powers of the 
inverse mass, and N. Therefore even in this simple case, there could be an aging regime 
at intermediate times. If this is the case it is likely that the methods of the present paper, 
using the GVA, will provide a basis for describing this out of equilibrium regime. 

In the limit of a zero mass, this model corresponds to the problem of diffusion in a random 
potential which has been extensively studied. In the particular case 7 = N = 1 model (1.1- 
(1.5) is the celebrated Sinai's model 26 . More generally it is known that <f>(t) ~ log(t)( 2/, ( 1-7 ^ 
for any finite N . However the two time correlations and thus the aging properties have not 
been studied previously analytically, except in the case of an applied force in Ref Fevi (see 
also Ref Mapa). 

To summarise, in this paper we study analytically the long time behaviour of model 
(1.1) with general V in the simplest case of zero internal dimensions. We use the formalism 
of Refs. [10,11] and we concentrate on the low temperature phases. We also reproduce and 
explain as explicitly as possible, some of the necessary calculations. The paper is organised 
as follows. In Section II we present the general mean-field dynamical equations for the 
zero internal dimension case. In Section III we describe the separation of the asymptotic 
dynamics in a stationary regime for long but similar times and a non-stationary regime 
for long and very different times. We derive the dynamical equations for both regimes. In 
Section IV we review the extensions of the equilibrium dynamics theorems to the asymptotic 
out of equilibrium relaxation. In Section V we present the - time-reparametrisation invariant 
- equations for the aging regime. We also discuss the strategy we follow in Sections VI and 
VII to find their solutions. In Section VI we use a 'one blob' ansatz to solve the short-range 
models and in Section VII we use an ultrametric ansatz to solve the long-range models. 
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Finally, in Section VIII we present our conclusions. 
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II. MEAN FIELD EQUATIONS IN THE LARGE TIME LIMIT 



The general dynamical equations for the two-time response function R(t,t') and correla- 
tion function C(t, t') for the model defined by Eq. (1.1) in the limit N — > oo were derived in 
Ref. [12]. For the special case of the p-spin spherical model they were presented in Ref. [10]. 
The dynamical equations assuming time-translation invariance (TTI) at all times for all the 
two-time functions were previously presented in Ref s. [24], [25], [30] and [31] for model (1.1) 
with short and long correlation, the p-spin spherical spin-glass and the RSGM, respectively. 
In these references the equilibrium dynamics a la Sompolinsky has been studied. 

In Appendix A we sketch a derivation of the general dynamical equations using the 
Martin-Siggia-Rose formalism complemented by a Gaussian approximation that becomes 
exact in the N — > oo limit. The mean-field dynamical equations in terms of C and R 
defined in Eq.(1.4) read 

dR (y) = R( t , f) + 4 f ds V"(B(t, s)) R(t, s) (R(t, f) - R(s, t')) (2.1) 
at Jo 

dc{ t tr) = -» + 2 f ds v ' W> s )) s ) 

at Jo 

+4 /* ds V"{B(t, s)) R(t, s) (C{t, t') - C(s, t')) + 2T R(t', t) . (2.2) 
Jo 

We use the Ito convention lim e ^ R(t,t — e) = 1, R(t,t) = 0. 

Using the fact that the last term in Eq.(2.2) vanishes for t > t' one also has 
ldC(t,t) 



-fiC(t, t) + T + 2 /* ds V'(B(t, s)) R(t, s) + 
Jo 



2 dt jo 

2 /* ds V"{B(t, s)) R(t, s) (C{t, t) - C(s, s) + B(t, s)) . (2.3) 
Jo 

For these problems it is convenient to use the mean squared displacement correlation 



function B(t,t') = C(t,t) + C(t',t') - 2C(t,t') = ((<f>(t) - ^{t')) 2 ) that satisfies for t > t' 
laB(M') V 



m 2 [C(t, t) - C(t', t') + B(t, t')] + 2j o ds V'(B(t, s)) (R(t, s) - R(t', s)) 

+T + 2 /* ds V"{B(t, s)) R(t, s) (B(t, s) + B(t, f) - B(s, t')) . (2.4) 
Jo 
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Note that one can use B(t,t') (t > t') rather than C(t,t') since the Eqs. (2.1), (2. 3) and 
(2.4) are a complete description of the problem together with the identities B(t,t) = 0, Vt. 
In these equations there are two external parameters T and //. 

In deriving these equations a mean over random initial condition cf>it = 0) (with a Gaus- 
sian distribution of variance C(0,0)) is implicit. One could certainly obtain the dynamical 
equations for a particular initial condition cf>it = 0) = cf> or a different distribution of initial 
conditions and attempt to study their effects in detail. We shall not do so here but we shall 
restrict our analysis to the study of the above dynamical equations. 

The time-dependent energy-density 12 is 



In order to obtain spherically constrained models one has to, on the one hand, let // 
be a function of time fi(t) and, on the other hand, determine fi(t) imposing explicitly the 
spherical constraint. If C(t,t) is set to q, then Eq.(2.3) implies 



In particular, if q = 1 and V is given by Eq.(1.7) and we compare with the equation 
determining the Lagrange multiplier zit) that enforces the spherical constraint in the p-spin 
spherical spin-glass 10 , we obtain 



The dynamics is described by the set (2.1)-(2.4) of non-linear coupled integro-differential 
equations that admit a unique solution for each 'initial auto-correlation' C(0,0). The dy- 
namical set of equations is causal. One can then use a numerical algorithm to iterate the 
equations and construct the solution step by step in time. A numerical analysis of this type 
of the model (1. !)-(!. 5) with 7 < 1 has been carried out in Ref.[12]. 
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III. STATIONARY AND N ON- STATIONARY DYNAMICS 



In the high-temperature phase both time translational invariance (TTI) and the fluc- 
tuation dissipation theorem (FDT) hold. Eqs. (2.1)-(2.4) reduce to a single equation. In 
Refs. [24,25] the model (1.1)-(1.5) with short and long-range correlation have been studied, 
respectively. The detailed analysis of the resulting equation shows that there is a critical 
curve T c (fi) below which no solution with these characteristics exists. This marks the end 
of the high-temperature phase. Similar situations occur in the p-spin spherical model 30 and 
in the Random Sine Gordon model 31 . 

The analytical study of the mean-field dynamical equations in the low-temperature phase 
requires the use of certain hypotheses. For a purely relaxational dynamics like a Langevin 
process that approaches an equilibrium situation one can show that the auto-correlation 
function is a monotonous function of the time difference r and that the response function is 
a positive function of its argument r. For a general out of equilibrium process the property 
of monotonicity is not obvious. However, for the kind of models we study it is a quite natural 
starting point to assume that the systems continuously drift away. As a consequence, though 
the auto-correlation function is not an exclusive function of the time-difference r for all times, 
it is assumed to increase monotonously when the separation between the two times increases: 



The weak-ergodicity breaking hypothesis 20 ' 10 includes the assumption above and, if the par- 
ticle moves in an infinite dimensional space and there is a finite mass: 






linw,£(M') 



k 



O 



V f 



finite , 




lim T 



oo 



lim, 



oo 



T + *',*') 



= q 



and 
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If the mass is zero from the start the particle diffuses and, in principle, both b and q are 
non-trivial time-dependent functions of t that tend to infinity when t — > oo. 

Hence, in the large-times limit, one can separate the relaxation that takes place for 
large and similar times such that it — t') jt' <C 1 from the large and widely separated times 
relaxation such that it — t')/t' = 0(1). The nature of the relaxations in these regimes is 
clearly different. The former is given by an 'equilibrium' or stationary dynamics while the 
latter is given by a 'non-equilibrium' or non-stationary dynamics with its maybe associated 
aging effects 10-12 . This hypothesis has been verified for different models with simulations, 
numerical analysis and exact results. For instance, the Monte Carlo simulations of various 
finite dimensional and mean-field spin-glass models 32 , the numerical study of the model 
defined by (1.1) for long-range correlations 12 , and the exact analytical results 13 for the 
model defined by Eqs. (1.6)-(1.7) when p = 2 provide evidence for this scenario. 

In the large t' limit we thus separate the t dependence in two regimes: for small time- 
differences it — t')/t' <C 1 there is equilibrium dynamics obeying the FDT and TTI: 

lim R(t + t',t') = r F (r) , lim B(t + t' , t') = b F (r) , lim C(r + t' , t') = c f (t) . (III.4) 

t'—^co t'—^co t'—^co 

The response and displacement functions are related by the FDT theorem 

= h dJ wr " M = 4 °-ir * M • < IIL5 > 

For large time-differences it — t')/t' = 0(1) there is slow dynamics violating both the 
FDT and TTI: 

R{t, f) = r(t, f) , B(t, f) = b(t, f) , C(t, f) = c(t, f) . (III.6) 

The limiting values for b are then 

b F (t,t) = lim^oo b F (r) = q 

(111.7) 

b(t,t_) = q b(t,0+) = b 

(see Eqs.(III.2)-(III.3)). We denote the initial time of the asymptotic regime s = + . It has 
to be understood as a very large time, as opposed to s = which is the actual initial time 

15 



and hence finite. One thus has b(t,0) > b(t,0 + ) = b . For fixed t the aging regime extends 
from t' = 0_|_ to t' = t- and FDT regime from t' = t_ to t' = t. We shall use lim^oo t-/t = 1. 
Calling q 0} qi and q the large time limiting values of C 

C(t,t) = c F (0) = q lim T ^ 00 c F (r) = qi , 

(III. 8) 

c(t,t_) = q x , c(t,0+) = q , 

and it follows from these definitions that 

q = 2(q- qi ) b = 2(q-q ). (111.9) 

Note that the corresponding definitions in the statics are as follows 5 ' 17 . One defines 
lab = ~N(4 > a4 > b) wnere a -, b are replica indices. The equal time correlation q stat corresponds 
to q aa . Within a replica symmetry breaking solution q a b for a ^ b is parametrized by q(u), 
< u < 1. While the precise shape of q(u) depends on the model, one has generally a 
constant q(u) = qf at from u c < u < 1 and q(u) = qQ tat . 

The second and crucial assumption can be called weak long-term memory 10 . It states 
that the response to any small perturbation applied during a finite interval [0,1^] eventually 
decays to zero for a large enough subsequent time t. In other words, the integral of the 
response function R(t,t') over t' in interval t' £ [^1,^2] with t 2 — t\ finite vanishes for long 
enough t. In a spin system this means that the thermo-remanent magnetisation decays to 
zero for a long enough time after having switched-off the (small) magnetic field. 

This separation in stationary and non-stationary dynamics allows us to write the dy- 
namical equations (2.1)-(2.4) for the two regimes, separating cleanly each contribution. In 
what follows we shall refer to the former time-regime ((t — t')/t' <^ 1) as the FDT-regime 
and to the latter time-regime ((t — t')/t' = 0(1)) as the aging regime. 

In Section V we shall present some numerical evidence for the separation in these two 
distinct time-regimes. 

The FDT regime 
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The equations obeyed by the FDT part ((t — t')/t' <C 1) were obtained in Ref. [12]. We 
recall and discuss them here lor completeness. 

dbp(r) 



dr 



drp(r) 
dr 

c F (0) 



2T + (-// + M) b F (r) 

-4 [ T dr' V"(b F (T - r')) r F (r') 6 F (r') 
Jo 

-4 / dr' [ (V"(b F (r + r')) r F (r + r') - V"(6 F (r')) r F (r')) b F (r') 
Jo 

+ (V'(6 F (r + r')) - V'{b F {r'))) rp(r') ] , (111.10) 

(-// + M) r F (r) - 4 T Jr' V"(6 F (r - r')) r F (r - r') r F (r') , (III.ll) 

Jo 

1 



Wl 



ith 



\i - M 

T + 2 / dr 1 V"(6 F (r')) r F (r') 6 F (r') + 2 / Jr' V'(6 F (r')) r F (r') , (111.12) 

v JO JO / 



M = 4 lim T V"(5(*, s)) s) . (111.13) 

t^co J 

M F = 4 / Jr' V"(6 F (r')) r F (r') , (111.14) 
Jo 

~M = M — M F = 4 lim / V"(6(t, s)) r(t, s) , (111.15) 

t^oo Jo 



It is easy to prove that thanks to the FDT relation (III. 5) the first two equations collapse 
and yield only one equation involving bp. However, it also contains contributions from the 
aging regime through a single, a priori unknown quantity 'the anomaly' M. We shall see that 
conversely the equations written in the aging regime do not include the details of the FDT 
solution. Thus, in order to find the full solution for the FDT regime the correct procedure 
is to first obtain the solution for the aging regime, compute the anomaly, insert it back in 
the FDT equation, and then solve for the FDT regime. 

The equation for the FDT regime is formally identical to the equation one finds studying 
the equilibrium dynamics a la Sompolinsky though the meaning of the anomaly is different 
in the two contexts. In the out of equilibrium situation the anomaly represents the memory 
of the system and, as can be clearly seen from its definition (III. 15), it is associated to an 
integration from the initial time of the asymptotic time-regime up to the final time t — > oo. 
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In the equilibrium dynamics instead, the initial time is chosen to be — oo in such a way to 
propose that the system is at some equilibrium state at a finite time and that it then is 
able to visit different ergodic components at diverging (with N) time scales. The anomalous 
term is then the contribution from all times starting at — oo and it is related to contributions 
from barrier crossing. 

Despite the different interpretation for the anomaly one can still borrow or at least 
formally compare some of the results concerning this equation from the previous works in 
which the equilibrium dynamics has been studied. The equilibrium analysis has been done 
in great detail in Refs. [24,25] for the short and long-range model (1.1) and in Ref.[31] for the 
RSGM. In the high temperature phase one finds M = 0. However, for T < T c (fi) the FDT 
equation has no solution for arbitrary large time differences if one keeps M = 0. There is 
a critical time-difference scale beyond which the FDT solution becomes unstable. It is then 
assumed that the equilibrium dynamics is 'marginal' meaning that the anomaly is chosen 
to have the minimal value required to remove the instability of the FDT solution. 21 ' 24 ' 25 

Here instead we shall first solve for the well-defined aging solution, compute the anomaly 
and then return to the study of the asymptotic (large time-difference r, i.e. lim T _ ) . 00 lim T / t /^ ) 
of the FDT regime. It turns out that, for the models studied here, the value of M obtained 
by this method coincides with the one obtained from the postulate of 'marginal stability'. 
Thus the detailed behaviour in the FDT regime, such as the power law decay of displacement 
to the asymptotic value q also coincides. 

The aging regime 

In Appendix B we derive the set of equations for large and widely separate times, t > 
t' 3> 1 and it — t')/t' = 0(1). The self-consistent procedure is as follows. In this slowly 
varying time-region, the time-derivatives in the l.h.s. of Eqs.(2.1)-(2.4) are assumed to be 
small compared to the r.h.s. in the equations and they are then neglected. This is related 
to the weak ergodicity breaking hypothesis that states that, as time elapses, the relaxation 
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of the systems is slowed down 20 ' 10 . Typically, the integrals in the r.h.s. go from s = to 
s = t or s = t' . In order to approximate them we start by separating the contribution 
of finite times, namely integrals going from the actual initial time s = to the starting 
time of the large time regime s = 0+. These are assumed to be small and hence neglected. 
This assumption is related to the hypothesis of weak long-term memory 10 and it allows us 
to write the long-time equations exclusively in terms of the long-time functions b and r. 
Finally we separate the large time-integration intervals to distinguish the FDT and aging 
contributions. We are typically left with integrals of products of slowly and rapidly varying 
functions. We approximate then the integrals by replacing the slowly varying functions over 
the whole integration interval is just constant an equal to the value at the border. 
The large-times dynamical equation for the response function r reads 

= r(t, f) (-p + 4 J* ds V"(b(t, s)) r(t, s) - | V"(b(t, *'))) 

-4 /* ds V"(b(t, s)) r(t, s) r(s, f) (111.16) 
■h' 

and that for the displacement correlation function b 

= b(t, t') + 2 J* ds V"(b(t, s)) r(t, 5)) + T + 1(V'(?) - V'(b(t, 0) 

+2 /* ds V'(b(t, s)) r(t, s) - 2 /* ds V'(b(t, s)) r(t', s) 
Jo Jo 

+2 /* ds V"(b(t, s)) r(t, s) b(t, s) - 2 /* ds V"(b(t, s)) r(t, s) b(t', s) 
Jo Jo 

-2 /* ds V"{b(t, s)) r(t, s) b(s, f) . (111.17) 
Jt' 

(NB. We have dropped all the subindices — , + in the limits of the integrals. In particular, 
we have dropped the + subindex in the zero limits. The zeroes and all the (time) integration 
variables should be interpreted as being in the asymptotic time regimes.) We shall call these 
two equations the r-Eq. and the 6-Eq., respectively. 

The equation for C(t,t) gives, at large t, another equation which contains contributions 
from the slowly varying parts as well as the FDT parts. In this large times limit, C(t,t) is 
assumed to have reached its asymptotic value q. Thus, 
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= - fl q + T + ^V(q) + 2 [ ds V'(b(t, s)) r(t, s) + 2 [ ds V"(b(t, s)) r(t, s) b(t, s) . 
1 Jo Jo 

(111.18) 

For completeness we give the equation lor c(t,t'). Using similar methods one has: 

= c(t, f) (-n + 4 J* ds V"(b(t, s)) r(t, a )) + |V W, *')) 

+2 / ds V'(b(t, s j) r(t', s) - 4 / ds V"(b(t, s j) r(t, s)c{t\ s) 
Jo Jo 

-4 /* ds V"(b(t, s)) r(t, s)c(s, f) . (111.19) 
Jt' 

We shall call it the c-Eq. 

We can obtain some relations between the asymptotic values q } q } b by considering par- 
ticular values of the times t, t'. Letting t' — > t_ in the above equations and using the limiting 
values (III.7)-(III.9) we find two conditions. Indeed, assuming that the integral in the r-Eq. 
vanishes in this limit, the condition for r(t,t_) to be non-zero is that the bracket in the 
following equation vanishes: 

= r(t, f _) (-fji + 4 J* ds V"(b(t, s)) r(t, s) - ^ V"(?)) . (111.20) 

The solution r(t,t_) = corresponds to the high-temperature phase for which there is no 
out of equilibrium dynamics. The equation arising from the requirement of the vanishing 
bracket involves the anomaly: M is equal to the integral above when t — > oo (see Eq.(III.15)). 
Hence, in order to have a non-trivial low-temperature solution, the anomaly must satisfy 

M = p + ^V"(q). (111.21) 

This equation is also the marginality condition that determines the anomaly in the equi- 
librium dynamics treatment. We see here how it arises naturally in the off-equilibirum 
approach. In addition it is related to the condition of vanishing replicon-eigenvalue in the 
static - replica - approach (see Section VI). 
Similarly, the 6-Eq. when t' — > t_ yields 
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2T r f 

= -// + — + 4/ dsV"(b(t,s))r(t,s) (111.22) 
^ Jo 

Subtracting Eq. (111.22) from Eq. (111.20) yields 

q 2 V"{q) = -T 2 . (111.23) 

This equation determines q as a function of T and the potential correlation V. Interestingly 
enough the equation is valid for all the potentials for T < T c . It is also independent of 
the mass and of the dynamical solution of the model and hence q is a purely geometrical 
quantity related only to the potential correlation. Besides, it imposes a condition on V": it 
must be negative at q in order to let (III. 23) have a sensible solution for q. The remaining 
equation, Eq. (III. 20), is important to select the behaviour of the models in the aging regime. 
Combining both equations one obtains: 

2T 

'/ = =fj (IIL24) 

fi — M 

One does not obtain any new equation by letting t' — > t_ in the c-Eq. One can check that 
when combined with Eq. (III. 18) it gives back Eq. (III. 22), as expected. 
Letting t' — > in the 6-eq. one finds: 

° = (" 2 + 2 Jo dS V ' W ' S)) S) ) K + T+ T {V{q) ~ V{ho)) 

+2 /* ds (V'(b(t, s)) + V"(b(t, s)) b(t, s)) r(t, s) - 2 f ds V"{b(t, s)) r(t, s) b(s, 0) , (111.25) 
Jo Jo 

b = b(t,0). This is the equation that fixes b . It is important to note that this equation 
exists only for a strictly non-zero mass. 

Finally, the energy-density (2.5) can be expressed using the separation of FDT and aging 
regimes as 

£(t) = \ q + | (V(0) - V(q)) - 2 J* ds V'(b(t, s)) r(t, s) . (111.26) 
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IV. EXTENSIONS OF THE DYNAMICAL THEOREMS TO THE NON-EQUILIBRIUM 

REGIME 



In this Section we review the extensions of the equilibrium theorems for the out of 
equilibrium dynamics of systems with slow relaxations proposed in Refs.[10,ll]. Since we 
shall mainly use in our calculations the displacement instead of the correlation function we 
here reexpress these extensions in terms of the displacement function. As opposed to the 
case of mean-field spin-glass models, in this paper we deal with problems that do not have 
a normalised correlation function. However, if the mass is non zero, the quadratic potential 
associated to it ensures that the equal-times correlation function reaches a finite limit at 
large times. This indeed can be easily checked solving numerically the dynamical equations. 
Hence, in these cases there is basically no difference between working with the displacement 
or with the correlation. We shall not discuss here in detail other further extensions associated 
to the massless cases in which there is unbounded diffusion and the equal-times correlation 
functions do not reach asymptotically a finite limit 16 . We finally discuss the description of 
possible singularities at the extremities of the aging regime. 

Definitions 

One defines the function X(t,t') as: 



with t > t'. In the FDT regime X(t,t') = Xp = — 1/(2T) while in the aging regime it 
measures the deviation from the FDT theorem. 

The mean squared separation B(t,t') monotonically increases when the two times t and 
t' become more separated. It increases from B(t,t) = to its maximum value B(t,0), that 
explicitly depends on the initial value for the correlation C(0, 0): B(t,0) = C(t, t) + C(0, 0) — 
2C(t, 0). Since we shall only consider large times, the maximum value that Bit, t') can take 
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is then B(t,0 + ) = b(t,0 + ) = b . Note however that, in general, B(t,0) is different from 
6(t,0 + ), more precisely, B(t,0) > b(t,0 + ). The monotonicity of B(t,t') allows us to invert 
the function B(t,t') and to write the function X, for instance, in terms of the smaller time 
and the displacement function, X(t,t') = X(B(t,t'),t'). One then assumes that in the large 
t' limit X(B(t,t'),t') approaches a function that depends only on B(t,t'), i.e. 

X(t,t') ~ X(B(t,t')) for t' > f . (IV.2) 

We shall also assume that is a monotonously decreasing function of B. In Section 

V we shall present results from the numerical solution of the dynamical equations that 
supports the assumption made in Eq.(IV.2) and the further assumption of monotonicity. 
The property of monotonicity permits to invert X(B) to have a function B(X). This 
function is useful to compare the asymptotic dynamical values with the equilibrium ones 
and it is somehow related 10 to the - replica - Parisi function b(u) = 2(q stat — q(u)). 

Similarly for three given times, thanks to the monotonicity of the displacement B, one 
can always write the relation 

B(t,t') = f(B(t,s),B(s,t'),t'). (IV.3) 

One then assumes that when all times are large, t > s > t' the explicit time dependence 
disappears and 

B(t, f) = f(B(t, s), B(s, t')) t > s > t' -> oo . (IV.4) 

The function / connects any three correlation functions; it has been called a triangle relation. 
Formally, one can also define the reciprocal function f(B' } B) such that 

B(s,f) =J(B(t,s), B(t,t')). (IV.5) 

Similar and in the end equivalent definitions of X and / could be given in terms of 

C(M0- 
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Properties of the triangle relations 



Let us recall the general properties of the function f(x, y). By definition f(x, y) is asso- 
ciative f(f(x,y),z) = f(x,f(y,z)). Since b(t,t') is an increasing function of the separation 
between t and t' one must have 

f(x,y) > max(x, y) . (IV. 6) 

One defines 'fixed points' b* such that /(&*, b*) = b*. Clearly f(x, x) > x for x between two 
fixed points. One can show that 

b* < x /(&*, x) = f(x, b*) = x , (IV.7) 
x < b* f(x, b*) = /(&*, x) = b* . (IV. 8) 

Hence the relation between a fixed point b* and any other point x is in this case 'maximum'. 
Between two fixed points b\ < x < b* 2l one can show, under the assumption that f(x } y) is 
smooth enough 34 , i.e. that it has a formal series expansion, and using the fact that there 
exists an identity (6J) and a zero (b^): 

i. f(x } y) = f(y } x) f is commutative 

ii. f(x,y) = J~ 1 (j(x)j(y)) f is isomorphic to the product 

with j(bl) = 1 and jib^) = 0. This implies y = f(x, z) = j -1 (j(z) / j(x)) } z > x. 

It is useful to explicitly compute some derivatives of f(x } y) and f(x,z) w.r.t. x, y and 
z when they are isomorphic to the product and division. For x, y £ (6J, b^) one finds 



df(x,y) 



dx 



df(x,y) 



x=y- 



dy 



1 . (IV.10) 



x=y- 



Moreover, since f(b* } x) = x and f(b* } x) = x for all x > b* then df(b* } x)/dx = 1 and 
df(b* } x)/dx = 1 for all x > b*. Instead, f(b* } x) = b* and f(b* } x) = b* for all x < b* } then 
df(b* } x)/dx = and df(b%x)/dx = for all x < b*. Thus, df(b*,x)/dx = 9(x - b*). 
It is useful to note that the form (IV.9-ii) can be expressed as: the form: 
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B(t,t')=j 




(IV.ll) 



where t\ is an arbitrary fixed time such that t\ < t' < t. 

However, the requirement that / has a formal series expansion 34 turns out to be too 
restrictive and for certain dynamical problems the large time dynamical equation may not 
admit a solution with an analytic j or j -1 in the whole interval b^). 

Indeed, studying the p = 2-spherical spin-glass model 13 one finds that the exact solution 
to the dynamical equations can be written, in the regime of large and widely separated 
times, as a function / isomorphic to the product as in Eq.(IV.9), but with a not well- 
behaved j(A) when A — > 1 . This is indeed related to the fact that the time derivatives 
dc(t,t')/dt = dc(t,t')/dt' = when t' — > t_ that implies that the inverse is not defined in 
this limit. 

More precisely, the zero temperature solution corresponds to 13 



with A = hit') I hit) = t' /t. j -1 has a vanishing derivative when A — > 1 and hence j'(A)| A ^ 1 — > 
oo. A — > 1 for j -1 corresponds to 6' — > 6 for / and / or t' — > t_ for c and b; i.e. it is the 
beginning of the aging regime. 

One of the main result of the present paper is that many of the models studied here do 
not admit a smooth solution when b' — > b but, instead, are solved in the aging regime by an 
ansatz still of the form (IV. 9) with a non-analytic j -1 at A — > 1. 

General organisation of fixed points 

We concentrate on the analysis of the long-time dynamics for which the functions defined 
above have been proposed. Having neglected the time-dependence in the function X and / 
suggests that the evolution of the system should be measured in terms of the displacement- 
value (or the correlation- value) instead of in terms of the times. This implies that time-scales 
are replaced by correlation-scales. A correlation scale is the range of correlations between 



j- 1 (X) = 2V2 



A 3 / 4 
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two fixed points of the function /. We call a 'blob' a discrete correlation scale between two 
separate fixed points, for which the function / is proposed to satisfy Eq. (IV. 9). 

It is clear that B(t,t) = is a fixed point of / for all the models we consider. The 
separation of large times into close and separated amounts to assuming the existence of a 
fixed point at Bit, tJ) = q that marks the end of the FDT correlation scale and the begining 
of the out of equilibrium correlation scales, the aging regime. 

When the times are close to each other, the displacements are homogeneous functions 
of time, B(t,t') = B(t - f). Calling n = t - t' , r 2 = t - t", r 3 = t" - t' and B x = £(ri), 
B 2 = B(t 2 ) } B 3 = B(t 3 ) one can immediately show that a function / as defined above exists 
in the FDT regime. Indeed, the displacement functions can be inverted to give r 2 = t 2 (B 2 ) } 
t 3 = t 3 (B 3 ) and then using n = r 2 + r 3 we have B x = B(t x ) = B(t 2 + r 3 ) = B(t 2 (B 2 ) + 
t 3 (B 3 )). Thus one finds that in the FDT regime B(t) = j _1 (exp(r)). As regards the 
function X it is just the constant associated to the FDT theorem in this correlation scale, 
X = -1/(2T). 

Another particular value of the displacement function is b . Since it is defined as the 
maximum value the displacement function can take in the large-times dynamics (for // > 
), it must be necessarily a fixed point of /. The displacement cannot go beyond this value. 
Moreover 0+) = 6(*',0+) = b . 

The strategy is to solve the large-times dynamical equations, viz. the r-Eq., 6-Eq. and 
c-Eq, (Eqs. (III. 16)-(III. 1 7) ) using the general properties of the functions X and /. The 
dynamical equations for each particular model, i.e. for each particular V(x), will determine 
the form of X and /. 
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V. REPARAMETRISATION INVARIANT EQUATIONS FOR THE AGING REGIME 



In this Section we write the dynamical equations for the aging regime, Eqs.(III.16)- 
(III. 19) in terms of the functions X and / defined and decribed in the previous section. We 
discuss the implications of those assumptions and we also present numerical evidence for 
them. We finally describe the strategy we shall follow in the two next Sections to find an 
analytic solution to the non-stationary dynamics. 

Neglecting the time-derivatives in the dynamical equations and the explicit dependence 
on time of the functions X and / amounts to introducing an artificial time-reparametrisation 
invariance into the large-times dynamical equations. The solutions we shall find are as 
a consequence invariant under time reparametrisations. This problem has already been 
encountered when solving the asymptotics of mean-field spin-glass models. The question on 
how to select the unique actual solution from the time-reparametrisation invariant family 
of solutions has not been answered yet. Indeed, this is already a well-know problem in 
the theory of non-linear one time differential equations. It is sometimes called the matching 
problem and it has been solved only for some particular cases. How to match the solutions we 
find for short time differences with one representative of the family of large time-differences 
solutions remains as a hard open problem. 

We can thus completely eliminate the times in the r-Eq., 6-Eq. and c-Eq, and rewrite 
them just in terms of the displacement b. With this aim we define 



The function M[b] is related to the anomaly: M = M[b ], see Eq. (III. 15). These functions 
could be also defined in such a way to contain the FDT regime. The upper limit in the 
integrals should then be set to b' = 0. In this case the second integral would be related to 





M, Eq.(III.13). 

We can then rewrite the r-Eq., Eq.(III.16), as 
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= -Hp + M) F[b(t, t')] + X M[b(t, t')] 

-4V"(6(*, 0) K*, F[q] - 4^ f V"(6') &(*, *'))] . (V.2) 

at Jb(t,t') 

Using F[q] = and integrating over t' one finds: 

= F[b] + M) + X M [6 ] + 4 £ d6 ' V "(6') X[6'] F[J(b', b)] (V.3) 

the integration constant being zero as can be seen by evaluating the resulting equation at 
b = q. This equation is in fact equivalent to a first integral of the r-Eq. 
The 6-Eq., Eq.(III.17), can be rewritten as: 

=(-// + M) b + 2T + ^ (V(q) - V(b)) + 4V'(6o) F[b' ] + 4 f db' V'(fc') X[b'] + 



■J 00 
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4 / db' V"(b') F[f(b, b')] + 4 / db' V"(b') b 1 X[b'] - 4 / db' V"(b') f(b, b') X[b'] 

J bo J bo J bo 

+4 [ b db'V"(b')J(b',b)X[b'], (V.4) 



where we denote b' = b(t\ 0). For t > t' > 0, bit, 0) > , 0). Since in the large time t' limit 
b(t', 0) = b' must be a fixed point, for all t' > 0: b = b(t, 0) = f(b(t, t'), b(t', 0)) = b(t', 0) = 

Both equations (V.3) and (V.4) are evidently time-reparametrisation invariant. In the 
same way we can re- write the c-Eq., Eq.(III.19), in terms of b. 

The 6-Eq. (V.4) evaluated in b — > q and b — >■ b gives back Eqs. (III. 22) and (III. 25), 
respectively, that in terms of the functional F and M read 

= (-n + M)q + 2T , (V.5) 
= -fib + 2T + ^ (V'(?) - V(b )) + 4V'(6 )F[6 ] + 4 £ A (V'(fc') b') X[b'] , (V.6) 

again, the latter exists only if // ^ from the start. 

The solution to the aging regime amounts, in this language, to solving Eqs. (V.3) and 
(V.4). Our strategy is as follows: 

A. We propose the existence of only one discrete correlation scale (one blob) with a constant 
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X[b] = x. We propose a completely regular /; by this we mean that the integration of a 
finite function of b' times /(&', b) between b and q tends to zero when b — >■ q. We find that 
only the p-spin spherical spin-glass and the limit of vanishing mass of the RSGM are solved 
by an ansatz of this type. 

B. We propose the existence of only one discrete scale but we allow for non- analytic solutions 
around b ~ q. This implies, in particular, that the integral above contributes in the limit 
b — >■ q. We find the conditions the potential correlation must satisfy in order to allow for 
such a solution. Model (1.1) with short-range correlations and the massive RSGM are solved 
in this way. 

C. We propose an ultrametric ansatz for all the correlations in the aging regime, namely 
/(&', b) = max(6', 6), b > q. We show the condition the potential correlation must satisfy 
to admit this ansatz as a solution. Model (1.1) with long-range correlations belongs to the 
family of models with an ultrametric dynamical solution. We present this analysis in the 
Section VII. 

In principle it is posible to show that the dynamical equations do not admit another type 
of solution with, for example, two discrete scales (two blobs) and different constant x inside 
each of the blobs or many other possibilities. This is a rather painful task that we do not 
do in this paper. We just study the three cases we describe above. 

A. Numerical evidence for the two regimes 

The numerical solution of the full set of dynamical demonstrates the existence of the two 
separated time-regimes and also the existence of a function X(B) for large times. One con- 
venient way to show this is by plotting the equivalent of the thermo-remanent magnetisation 
of a magnetic system 
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m(M') = / dt"R(t,t") (V.A.I) 
Jo 

vs the displacement function B(t,t'). If R(t,t') = X(B(t,t')) dB(t,t')/dt' for large t', then 
X(B) is the local slope of these curves. 

In Fig. 1 and 2 we plot m(t,t') vs. B(t,t') for the power law model with long 7 = .5 
and short 7 = 1.5 correlations, respectively. The external paramters are // = 0.2, 6 = 2.5 
and T = 0.2 and we choose C(0,0) = 0. Three curves for three times t are included in 
both figures: t = 200,300,400 (1000, 1500 and 2000 iterations of the algorithm with a step 
h = 0.2; all our numerical results have been performed with h = 0.2). They show clearly 
two B regimes: for B < q the slope is constant and equal to — 1/(2T) = —0.25 (FDT holds). 
For B > q and if the potential correlation is long-range the slope is not constant while if it 
is short-range the slope is still compatible with a constant \x\ < | — 1/(2T)|. The breaking 
point q and end point b of the displacement can be read from the figures and as we shall 
see in Sections VI and VII they coincide nicely with the analytical predictions. 

The curves for different total times t do not exactly lie upon each other showing that 
the numerical X(B) still depends on t. We expect that for sufficiently large times they will 
converge to the analytical prediction that we describe in the Sections VI and VII and that 
we include in these figures as the broken lines. We can see that the numerical convergence is 
indeed rather slow. In the insets we plot the same curves for the same choice of parameters 
but for the high-temperature phases, T = 1. t = 400. In both figures the slope is —.5 for 
all B. 
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FIGURES 



E 




q=1 .03 B(t,t') bo=6.05 

FIG. 1. m(t,t') vs. B(t,t') for the long-range power law model. 7 = .5, fj, = 0.2, = 2.5. For 
the full lines / = 200, 300,400. In the main plot T = 0.2 while in the inset T = 1. > T c . Note that 
with the choice C(0,0) = one has automatically C(t,0) = and thus B(t,0) = C(t,t). 
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FIG. 2. m(t,t') vs. B(t,t') for the short-range power law model. 7 = 1.5, fj, = 0.2, 9 = 2.5. For 
the full lines t = 200,300,400. T = 0.2 while in the inset T = 1. > T c . q ~ 1.17, b ~ 3.17. 
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VI. SOLUTIONS TO THE DYNAMICAL EQUATIONS: SHORT-RANGE POTENTIALS 



In this Section we apply a one blob ansatz to the aging regime equations and we show 
that those associated to the short-range potentials are solved in this way. 

A. Asymptotic values: comparison with the statics 

We start by studying the asymptotic values q, b 0} and x. Though a deep understanding 
of the relation between the dynamical X[b] and the statical replica function x[q] is still 
missing it is by now clear that, at least formally, they play similar roles. We compare here 
the aymptotic values q and b , and the parameter x = X[b] inside the blob, to the static 17 
q sta \ b s o tat and x stat . 

Equations (III. 22) and (III. 25) with X[b] = x are identical to the ones obtained within the 
replica analysis 17 as the extremum conditions for the static q and b . This is a consequence 
of an algebraic relation between dynamical and static (replica) approaches 35 . Instead, the 
third equation, Eq.(III.21), dynamical in origin, is formally equal to the condition to have 
a vanishing replicon-eigenvalue of the replica treatment. It does not coincide with the third 
equation arising from the replica recipe, i.e. the one obtained maximising the free-energy 
w.r.t. x. As a consequence, the dynamical value for x does not coincide with the static one. 
This implies that also the values for q and b differ since x enters in their equations. The 
dynamical phase diagram turns out to be different from the static, the dynamical transition 
temperature is higher than the static critical temperature for a range of masses. Moreover, 
the dynamical asymptotic energy is higher than the equilibrium energy. This discrepancy 
has already been observed for the special case of the p-spin spherical spin-glass and it seems 
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to be related to a peculiar organisation of the metastable states of the TAP 37 free-energy 
landscape, with the existence of a threshold (higher than the equilibrium free-energy, below 
which infinite barriers appear 10 ' 38 ). A more detailed analysis along these lines is in order to 
determine if this is at the origin of the difference between static and dynamical asymptotic 
values in these models too. 

If the model is massive, the particle is constrained to move in a restricted region of the 
space and hence b and q are finite. If we take the limit //—*■() afterwards we shall see that 
b — > oo in such a way that fib remains finite. Instead the case of a strictly massless model 
is more delicate. Let us now obtain the equations that fix q, b and x for a massive model. 
Eq. (III. 23) determines q. When x is constant the equations that fix b and x follow from 
Eqs. (V.5) and (V.6) and they read 

V'(q) - V'(b ) fiT 



K- q 
b - q 



2q 



1 



'2T 



2Ta 



q 



while Eq.(VI.A.3) gives q: 



1 

q = - 



T+j; V(q) + 2x (V(q) - V(b j) + 2x(qV(q) - b V(b )) 



The anomaly M can be immediately computed: 



(VI.A.l) 
(VI.A.2) 



(VI.A.3) 



M = ix(V(q)-V(b )) 



(VI.A.4) 



and can be used to solve the FDT equations (see Section III). 
The asymptotic energy-density reads 

m = f q + y - v ^)) - 2x ( v ^) - v ^)) 



(VI.A.5) 



and thus it is determined by the asymptotic values q } b and x. Since x does not coincide 
with the static - replica - the asympotic dynamical energy differs from the static one. We 
shall show below some numerical results on this point (Fig. 7). 
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B. The fj, — T phase diagram 



In Fig. 3 we present the phase diagram for the power law model with 7 = 1.5 and 6 = 2.5. 
The equation for q, Eq.(III.23), determines a first critical temperature (independent of 
the mass) given by 



= ^-(q*) 2 V"(q*) 



(VI.B.l) 



with q* the value of q at which q 2 V"(q) attains its maximum, i.e. 

q*V'"(q*) 



2V"{q*) 



-1 



(VLB. 2) 



The glass phase cannot extend beyond T = Tj 1 ) since there is then no solution to Eq.(III.23). 



0.45 



T 




FIG. 3. Phase diagram for the power law model with 7 = 1.5 and 9 = 2.5. 



Just from their definitions b > q. This condition determines a critical line T( 2 \ji). 
Eq.(VI.A.l) when b — >■ q yields 

= ^ • (VI.B.3) 
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where q(T) is the solution to Eq.(III.23). The general equation for this line i 



is: 



'2T^ ■■ 2 



v "l^rJ = -T- (VLB - 4 » 

This same equation is found in the replica treatment as the condition to have a a vanishing 
replicon eigenvalue. Statically this determines the line below which the replica symmetric 
solution becomes unstable and hence the transition line from the RS solution to the 1RSB 
solution 17 . Eq.(VI.B.4) delimits a domain in the ji-T plane were the aging solution exists. 



When T — > we find a finite critical // = fi c = 2y |V"(0)| beyond which the above equation 

has no solution. When — > and the potential is short-range (i.e. V"(x) decreases faster 

than 1/x 2 at infinity) T^ 2 \ji) — > 0. 

For the power law model the explicit formula for Tj 2 )(//) is 

/ fry \ 1/(7+1) \ 

^) = f((^J -0) ■ (VLB.5) 

If 7 > 1 , T^(fi -> 0) -> and for 7 = 1 , T^ 2 \ji -> 0) = 1/^2. 

Returning to the general potential correlation and evaluating x from Eq.(VLA.2) at the 
critical line b = q + e, e — > 0, one obtains 

q V"'(q) 

x = x c = ^- — — 7-^- . (VI.B.6) 

with q = 2T^ I n on the critical line. One sees that as long as the function q 2 V"(q) is 
increasing (small values of q } q < q*) } one automatically has \x c \ < 1/(2T). Thus, even if 
the transition is continuous for q and b there is a jump in x, from the value x c below the 
critical line to the high-T phase value x = xp = — 1/(2T). 

The value x c = — 1/(2T) is exactly reached at the temperature = where q 

attains its maximal value. Equations (VI.B.l) and (VLB. 4) determine the tricritical point 
(/i*,T*), at which b = q and x = — 1/(2T) simultaneously. For //<//* the nature of the 
transition changes. The phase boundary is now determined by the condition x c = — 1/(2T), 
while 6 at the transition is no longer equal to q. Setting x = — 1/(2T) in the equations 
(VI. A. 1) and (VI. A. 2) one obtains b = 2T^ / ' \i at the transition and the following equations 
for the transition temperature: 
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V(q) - V = [ -^- - , (VI.B.7) 

V v- J q 2 

with q given, as always, by Eq.(III.23). The transition is now continuous in x but discontin- 
uous in q. rj 3 )(/i) is well above the transition line for the statics and it reaches a finite limit 
at fi — > 0, as opposed to T c — > when // — > in the statics. For the potential correlation 
(1.5) one finds that when // -> 0, (T^(fi)) 2 -> ^^i 2 ^ 1 ) 1 and 9 = 0/( 7 - 1). 

For model (1.1) with V given by Eq.(1.5) and 7 > 1, these conditions imply a phase 
diagram identical to the one described in Ref. [24] for the equilibrium dynamics. It is not 
clear why it should be a difference between static and equilibrium dynamic transitions (and 
furthermore between static and asymptotic equilibrium dynamic energies) especially if one 
remembers that the initial assumption of the latter approach is that the evolution starts at 
an equilibrium state. In the out of equilibrium evolution there is a priori no reason for the 
two transition temperature to coincide and indeed one knows this happens in real systems 
as glasses. We shall come back to this discussion later when we shall present the numerical 
results for the asymptotic energy-density vs. the static one. 



C. Equations within a blob with X[b] = x 

Using X[b] = x: 

F[b] = -x(q-b) M[b] = 4x (V(q) - V(b)) , (VI.C.l) 
and the first integral of the r-Eq., Eq. (V.3), and the 6-Eq., Eq. (V.4), read 

= (b-q) [-fi + ix (V'(q) - V'(b ))} + Aq (± + z) (V'(q) - V'(b)) 

+ 4x [ b db'V"(b')J(b',b) (VI.C.2) 



= b [-// + 4x (V'(q) - V'(b ))] + 2T + 4q^ + (V'(q) - V'(b)) 

+ Ax f db' V\b') 7(6', b) . (VI.C.3) 

Jq 
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We see that these equations coincide if we use // — M = 2T/q. We shall concentrate then 
on the r-Eq. (VI.C.2). 



D. Analytic ansatz at /' ~ / 



We can now show that a regular solution at t' — > t_ cannot exist for all the models we 
consider. In fact, if / is regular at b — > q one can evaluate Eq. (VI.C.2) and its variations 
with respect to b in this limit. If / and all its variations with respect to b are smooth 
when b — > q } one can set the integral in each equation to zero (the limits collapse and the 
integrands are regular by hypothesis) when evaluating at b — > q. However, if this is so, one 
gets new conditions on x,q and b that are not always compatible. 

More precisely, the first new equation arises from the second variation of the r-Eq. w.r.t. 
b. When b — > q it implies 

* = ^ , (VI.D.1) 

where we have used df(b' } b) / db\^i^^ q = 1 and we have set the remaining integral to zero. 
For general models, i.e. general Vs, this condition is not compatible with the x arising from 
Eq.(VI.A.2) above. 

In an alternative way, one can investigate the equation the potential correlation V has to 
satisfy in order to allow for a regular solution to the dynamical equation. Combining Eqs. 
(111.23), (VI.A.l) and (VI.A.2) with Eq.(VI.D.l) yields 

v 'M~nk) = V"M , (VLD/2) 

We can now consider the temperature T as a varying parameter and then q as the indepen- 
dent variable in a differential equation for V(q). Calling y(q) = V(q) and y = V'(b ) we 
have: 
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-(V ~ Vo) W - (K - q)y") = (K - q) y' 2 . (VI.D.3) 

The potential correlation has to satisfy this equation in order to allow for a regular solution. 

We have checked that from the potential correlations described in the Introduction only 
the one associated to the p-spin spherical model with p > 2 admits a solution with a well- 
behaved j -1 when A — > 1 (the case p = 2 is particular because the first derivative of j -1 
is zero). The x arising from the second variation of Eq. (VI. C. 2) given by Eq.(VI.D.l) is 
compatible with the x arising from Eq.(VI.A.2). All the higher variations are trivial and 
do not give further constraints. Conversely, if one perturbs the original p-spin Hamiltonian 
with an arbitrary small term associated to a different p in the manner of Ref. [36], one 
immediately sees that a regular solution for the dynamics is not more allowed. One can say 
that the p-spin spherical spin-glass is somehow a marginal model. 

The problem for a general V is in fact so constrained that it is not hard to believe that 
the potential correlation associated to the p-spin spherical spin-glass be the unique function 
having a regular solution. 



E. Non- analytic solution at /' ~ / 

Using fi — M = 2T / q the remaining equation (VI. C. 2) reads 
= -l(b-q)-2q(x + ^j (V\b)-V\q)) + 2xj\b'V'\b')j-'(^j . (VI.E.I) 
It is useful now to remember the original time-dependence inside the integral using 

MMO^'gg. (V..E.2, 



10. 



Thus, defining a new variable u and a new function B(u) as 

exp(-u) = h(t) b(t, t') = q + B(u - u') = f 1 (^j^J , (VI.E.3) 
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Eq.(VI.E.l) becomes 

= B(u)V"(q) - V'(q + B(u)) + V'(q) - V"(q) / dB(u') V"(q + B{u')) B(u - u 1 ) , 

1 Jo 

(VI.E.4) 

or, integrating by parts the last term, 



= -B{u)V"{q) + V\q + B{u)) - V\q 

d 



+ ^V"(?) 



V'(q)B(u) + — / di/ V'(? + B(u')) B(u - u') 
Ou Jo 



(VI.E.5) 



This equation represents the homogeneized version of Eq.(VI.E.l). It is important to note 
that it is not homogeneous in time and thus its solutions are not either. The homogeneity 
holds at the level of non linear functions of the correlation functions and not at the level of 
times. It has no free parameters, q and x are fixed by Eqs. (III. 23), (VI. A. 1) and (VI. A. 2). 
The invariance of the aging regime equations under time-reparametrisations t — > hit) implies 
the invariance of Eqs. (VI.E.4) and (VI.E.5) under dilatations of it: u — > rju. 

Equation (VI.E.5) can be solved exactly for some special choices of V such as the p- 
spin spherical model and the RSGM in the limit of vanishing mass. The solutions are 
B(u) = 2qps (1 — exp(— rju)) and B(u) = rju } respectively. We shall describe these cases in 
detail at the end of this Section. In what follows we shall instead analyse the small and 
large u behaviour of B(u) for the general potential correlation. 

We solve Eq. (VI.E.4) around u ~ 0, i.e. b ~ q } that corresponds to B(u) ~ 0, by first 
expanding in powers of B(u) and then proposing a series expansion of B(u) in terms of u a . 
a is a positive exponent that should be determined by the equation. An exponent a smaller 
than one is related to an irregularity at u ~ 0. 

The expansion around B(u) ~ yields 

V {n+1) (q) (u( y , n , 2xq V<%) 3 



o = E 



n=2 



'^-1 T>( U -U') 



(VI.E.6) 



V"(q)n\ v x JJ T (n-l)\ du J 

It is easy to see that the zero-th and first order equations in B(u) are immediately 
satisfied. The second order equation is however non-trivial: 
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= ^ (B(u) f + ^ V\q)^- [ du< B(u')B(u - u<) . 



(VI.E.7) 




FIG. 4. B(u) vs. v}l a with 1/a = a 7 = 1.5, = 2.5, T = 0.2 and fi = 0. From top to bottom 
the number of coefficients are 7, 9, 10, 8. 



We now proceed by proposing a formal series expansion for the function B(u) around 
u ~ in terms of u a 



B(u) = B(u a ) = 



(VI.E.8) 



m=0 



and then by solving Eq. (VI. E. 6) term by term in the double series expansion. More 
precisely, Eq. (VI.E.7) at first order in u a implies 







y V"(q) T v ^ r(l+2a) J 
and this equation fixes the function 

(r(i + a )) 2 

p(a) = r(l + 2a) 
and hence the exponent a. The limiting values of p are 



(b lU c 



(VI.E.9) 



(VI.E.10) 
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lim^oo p(a) = , 

p(l) = 1/2 

lim„^ /o(a) = 1 , • 

For a general V 



B(u) 




(VI.E.ll) 



1 2 3 4 5 

FIG. 5. B(u) vs. u x > a where 1/a = a; 7 = 1.5, 9 = 2.5, T = 0.2 and /x = 0.2. The curve 
pointing upwards has 30 coefficients while the curves pointing downwards have 50, 20, 10 from left 
to right. The straight line is B = b — q. 

The first coefficient bi is not fixed as a consequence of the reparametrisations invariance. 
In Appendix D we describe how we can obtain the coefficientes b m} m > 2, in terms of 
bi and then reconstruct the series for B(u). In Figs. 4 and 5 we show B(u) for model 
(1.1) with 7 = 1.5, 6 = 2.5 and T = 0.2. Fig. 4 corresponds to the massless case and 
Fig. 5 to a massive case with p = 0.2. The different curves in both figures correspond to 
having approximated the series for B(u) with a sum with different numbers of terms. We 
see that both in the massless and massive cases B(u) can be obtained in this simple way 
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up to reasonable values of u a . We have checked however that in the massive case a Pade 
approximation gives B(u — > oo) ~ b — q as expected. In the massless case we obtain the 
analytic asymptotic behaviour of B(u) in the next section and find that it is in good accord 
with the tendency of the curves in Fig. 4. 



F. Solution for b ~ b , large u behaviour for the massless models 



The analysis of Eq.(VFE.4) or (VFE.5) for large u is delicate and we shall not pursue it 
in all generality. We shall describe instead what happens in the limit of a vanishing mass 
for a general potential correlation. 

In the limit // -> Eqs.(VI.A.l), (VI.A.3) and (VI.A.2) imply 

V'(,) - V'(W ~ -f- l>.i~'^- T -^\ (VI.F.12) 

x Zqx 1 (q) 

Thus b — > oo, V'(b ) — > and V'(q) ~ —T/(2qx). Using this result Eq.(VI.E.4) simplifies: 

V(q) - V(q + B(u)) = -^M |_ T du > V '( 9 + B(u')) B(u - u') . (VI.F.13) 

V (q) uu Jo 

Let us now study the asymptotic behaviour of B(u) at large u. Since B(u) — > oo one 
must have 

J du'V , (q + B{u'))B{u-u') ~ ^7^y M (VI.F.I4) 

Defining I = J °° du'V'(q + B(u')) and assuming I to be a convergent integral immediately 
implies B(u) ~ rju with rjl = • One can check that this is the only possibility for 

short-range models. It implies 

&( M ') ~ 9 + In (VI.F.15) 

in the limit hit 1 ) /hit) « I for all short-range models. This form which turns out to be 
exact for the RSGM in the massless limit (see next section) for all it, implies that r(t,t') is 
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only a function of t' for widely separated times. We have checked this result numerically. In 
Fig. 6. we plot the response function R(t,t') vs. the smaller time t' for four times t in the 
massless case. The parameters are T = 0.2, 6 = 2.5. 




600 800 1000 1200 1400 1600 1800 2000 
t' 



FIG. 6. The response function R(t,t') vs. t 1 for / 
fj, = 0,T = 0.2,7 = 1.5 



ih, h = 0.2 and i = 100,200,300,400. 



In the present Section we are letting //—*■() and using all the equations including (VI. A. 1). 
Note that an identical equation (and solution) for B(u) would be obtained by setting // = 
first and then setting b = oo abandoning the equation (VI. A. 1). 

A very interesting quantity to calculate for the short-range models in particular in the 
massless limit is the asymptotic energy-density 

£^ -+ 0) = fi q + | (V(0) - V(q)) - 2xV(q) . (VI.F.16) 

As we shall see in the following Subsection, the asymptotic energy-density of the massive 
power law model with short-range correlations (1.5) is higher than the equilibrium energy- 
density obtained with the replica analysis 17 . Moreover, within the dynamic approach we find 
a non-trivial solution when //—*■() with a finite asymptotic energy-density while statically 
the solution does not exist in the massless limit since T c (fi — > 0) — > in the statics. This 
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difference between asymptotic and equilibrium values seems to be characteristic of models 
solved by f RSB ansatz statically, and it is related to the fact that the dynamic B(X) is not 
equal to static Parisi function b(u). 



G. Some particular models 



The power law model 



For model (1.1) p(ct) becomes 



p(a) 



xq 27 



(0 + q) 



(VI.H.l) 



In Fig. 2 we show the analytical prediction (broken straight lines) for this model with 
7 = 1.5, fi = 0.2, T = 0.2 and 6 = 2.5. The full lines are the numerical solution to 
the dynamical equations for t = 200,300,400. from Eqs.(III.23)-(VI.A.2) q = 1.17 and 
b = 3.17. The FDT decay and q are quite rapidly reached by the numerical solution. 
Instead, the convergence towards the aging theoretical curve and b is slower. Note that the 
end of the full lines is at b(t, 0) that is slightly larger than b = 3.17, as expected. 

When the mass p tends to zero the asymptotic values q and b 0} and x can be easy 
obtained. They are given by the following equations 






7 




(VI.H.3) 



When fi — >■ 0, b t 



o 



00 and q — > 00 in such a way that b p and q/j, remain finite. Thus, 



in the zero mass limit 
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p(a) ~ i±I . (VI.H.4) 

Reconsidering the limiting behaviour of p } cfr. Eqs. (VI.E.ll), we see that 

7 — > oo — )■ 1/2 a — s- 1 

7 -> 1 =>• /? -> 1 =>• a ^ (VI.H.5) 
7<1^/?>1 a /3 

For 7 < 1 (long-range potential) p > 1 but, from its definition (VI.E.10), cannot be 

greater than 1. Hence, we have demonstrated that the long-range potentials of type (1.5) 
do not have a one blob solution, at least in the massless limit, for any temperature. This 
demonstration can be extended to show that indeed the long-range model does not admit a 
one blob solution even when the mass is finite. 

At the zero temperature limit we can obtain the explicit dependence of q, b and x on 
the parameters 6 and 7: 

9 ~T,A^< 1 x ~ -.fce 2 ^ =0(1). 

V7 V 2 (VI.H.6) 

bo, ~ ^ *V = 0(1) qp ~ ^-)/ 2 = 0(1) 

Note that x remains finite at zero temperature. 

The asymptotic energy-density for the power law model is 

£00 = M + j; (A 1 " 7 -(0 + q) 1 - 1 ) ~ x ^— y ((6 + q) 1 ^ -(6 + b,) 1 ^ 

(VI.H.7) 

in the massless case it reads 

- 0) = ^-{0 + q) + ^ j K" 7 - (* + 9) 1 " 7 ) • (VLH.8) 

In Fig. 7 we plot the energy-density decay for the massless power law model with parameters 
7 = 1.5, 9 = 2.5 and T = 0.2. Equation (VI. H. 8) implies ~ —.137. In the inset we 
plot In (£(t) + .137) vs. Int. One can see that the energy-density approaches its asymptotic 
value with a power law. 

46 








-1.2 - 
-1.4 - 
-1.6 - 



50 100 150 200 250 300 350 400 

t 

FIG. 7. Energy- density decay: 8{t) vs. t for the short-range model with 7 = 1.5, 9 = 2.5, 
T = 0.2. In the inset In (S(t) + .137) vs. Int. 

Finally, at zero temperature we have 




(VI.H.9) 



The zero mass random Sine Gordon model 

In the limit of zero mass, Eqs. (VI. A. 1) and (VI. A. 2) can be easily solved. The maximum 
displacement b tends to infinity in this limit but the product fib stays finite; we obtain 

q 2 exp(- 9 /2) = — - fib qn , (VT.H.10) 

A x q q 

Therefore, p(a) satisfies 

p(a) = ^ a = l. (VI.H.ll) 

Besides, this is a specially simple problem that can be solved completely from the original 
equations. Indeed, the solution is 
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b(t,t') = q + In 



( 



h(t') 



h(t) 



) 



r(M') 



4T h(t') 



q h'(t') 



(VI.H.12) 



The exponent a = 1 is consistent with this solution since it translates into 



B(u) = rju . 



(VI.H.13) 



The triangular relation function is thus f(x, y) = x + y — q. Note that a particular solution 
corresponds to h(t) = t s which we speculate is a good candidate for being the solution 
selected by matching at small times. 

The p-spin spherical model 

It is interesting to recall the known results for this model and compare them with those 
obtained with the general treatment described above. For spherically constrained models 
b = 2. In the case of the p-spin spherical model, with p > 2, Eq. (III. 23) implies p(p — 

1) /2 (1 — qps) 2 q P p~s = T 2 . The mass // in Eqs. (VI. A. 1) and (VI. A. 2) has to be interpreted 
as the constant long-time limit of the time-dependent fi(t) related to the Lagrange multiplier 
enforcing the spherical constraint, cfr. Eq. (2.7). Solving for x one gets x = — 1/(2T) (p — 

2) (1 — qps) I qps These equations coincide with those found in Ref.[10] Eq.(VI.E.4) is solved, 
for all it, by B(u) = 2gps(l — exp(— rju)) with rj undetermined. Hence the function j defined 
in Eq.(IV.9) that relates any three displacements in the aging regime is jib) = qps(l — b/2) 
and the triangle relation is just f(x } y) = xy/qps- 

If we study Eq.(VI.E.5) for u ~ as described in the previous section, from Eq.(VI.E.ll) 



The case p = 2 is particular. Both the third derivative in the numerator of Eq. (VI.E.ll) 
and x are zero. One cannot naively use Eq. (VI.E.ll) to determine a. Equation (VI. E. 5) 
can be used also to check that the large u analysis of it is delicate. 





a = 1 for p > 2 . 



(VI.H.14) 
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H. From short-range models to long-range models: The appearance of the ultrametric solution 

Let us study in detail the limit 7 — > 1 from above for the power law model. In general, 
we have solved the short-range models with an ansatz consisting of a constant x and a 
triangular relation f(x,y) = J~ 1 (j(x) j(y)), with J~ 1 (x) a non-analytic function at x = 1: 

j- 1 (I + e) = q + ce a + 0{e 2a ) , (VI.I.l) 

that is different for different models, i.e. different V. This implies for / 

f(q + Z} q + z')oc (z 1 '* + z a l a ) a + 0(z\ z'\ zz 1 ) . (VI.I.2) 

In Section VTH we have also shown when 7 tends to one the exponent a for the power 
law model tends to 0. Hence, when 7 — > 1 

f(q + z, q + z 1 ) ~ q + max(z, z 1 ) + h.o.t (VI. 1. 3) 

and we see explicitly how the ultrametric solution is approached when the potential corre- 
lation range is increased. 

The FDT theorem violation x stays constant when 7 — > 1 from above and this is also 
valid when the limit is approached from below, as can be seen from Eq.(VII.ll). 
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VII. SOLUTIONS TO THE DYNAMICAL EQUATIONS: LONG-RANGE POTENTIALS 



We shall apply the ultrametric ansatz to Eqs.(V.3) and (V.4). The ultrametric ansatz 



means 



b = f{b',b") = max(6',6") 
b" = J(b',b) = max(6',6) 

V 6, b" inside the interval [<?,& ]. Note that df(b',b)/db = 1 and df(b',b)/db' = 0, since 
b > V. Furthermore d^J \b' , b) / db^ = d (k ">J{b> \b) / ' db' (k) = 0, Vk > 1. In the dynamic 
equations (V.3) and (V.4) only / enters. We can then take variations of these equations 
w.r.t. b safely. 

Inserting the ultrametric form of / in Eq.(V.3) and in its variation with respect to 6, 
and using Eq.(V.5) yields 



= X[b] 



2q 
T 



(V"(q) - V"(b)) - 4M[6] + iV"(b)F[b] 



(vn.i; 



Solving these equations for F[b] and M[b] we obtain 



F[b] = T 



2T 



V"(b) - JV"(b)V"(q) 



M[b] = ±| 



V"{b) 

V"(q) - ^JV"(b)V"(q) 



(VII.2) 



The second variation of Eq. (V.3) combined with the above expresions for F[b] and M[b] 
gives 



q V"'(b) 

X[b] = 



V"(q) 



(VII.3) 



"4T V"{b) ^ V"{b) ' 

The sign must be chosen in such a way that X[b] is negative; it depends on the relative signs 
of V'" and V". In the following we shall consider V" < and V'" >. This is the case for the 
power law model with long-range correlations. We shall then choose the lower signs in the 
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expressions above. If / is a maximum, all the integrals appearing in the higher variations of 
the equation w.r.t. b in the limit b — >■ q vanish. 

We require that \X\ be a decreasing or constant function of its argument b. This indeed 
selects the models that can be solved by an ultrametric ansatz. Imposing |X[&]| < for 
b > q yields 



)/// 



V(b) _ V 



1 >Z7^ with v= , „, A3/2 • (VH.4) 



For instance, if the function V is of the form (1.5), we have 

fe + q ^ 



2 



< 1 . (VII.5) 



\e + bj 

If 7 > 1 this inequality cannot be satisfied. This means that model (1.1) with short-range 
correlations cannot be solved by an ultrametric ansatz. Instead, if 7 < 1, the condition is 
satisfied for all b and model (1.1) with long-range correlations can be solved by an ultramtric 
ansatz. 

As a second example we can examine the RSGM model. Eqs. (1.9) and (VII. 4) imply 

exp(6/4) ^ ^ ("VII 
exp(^/4) ~~ 

Since b > q then this imposes q = b and the RSGM does not admit an ultrametric solution. 

Finally, we study the p-spin model. In this case, l — B/2 = C and the correlation function 
decays from q when t' — > t to when t' — > 0. Then, 

< 1 (VII.7) 

Since q > c the inequality is not satisfied. 

In conclusion, from the particular cases we presented in the Introduction, only model (1.1) 
with long-range interactions, i.e. 7 < 1, admits an ultrametric solution. More generally, 
Eq.(VII.4) selects the potentials that admit an ultrametric solution. 



The formula (VII. 3) for X together with the ultrametric form of / represent the full 
solution for the model. Due to the time-reparametrisation invariance we have included in 
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the equations we cannot go beyond these expressions and obtain the explicit time-dependence 
of the displacement b. We can check that Eq.(III.17) is also solved by this ansatz. 

We can now explicitly obtain q, b and q and compare them with the static results 
obtained with the replica approach 5 ' 17 . Equation (III. 23) determines q and with X[b] given 
by Eq. (VII. 3) after a bit of algebra we obtain b from Eq. (V.5): 

V"(b ) = . (VII.8) 

For the power law model with 7 < 1 q and b are given by 

a 2 9T 2 11 

q (9 + 6 )"^ +1 )/ 2 = . (VII.9) 



and X[b] i 



is 



x[b\ = - l - (0 + b)W 2 . (vii.io) 



In the limit 7 — > 1 we have 

X[b] = * = (v11.1i) 

that coincides with the constant x obtained for the short-range models when 7 — > 1 + 

One can check that these equations for q and b coincide with the equilibrium equations 
obtained with a full replica symmetry breaking ansatz 17 ' 12 . One can also check that the 
expression we obtained for X(B) is formally identical to the Parisi function x(q) of the 
replica treatment. Hence, the asymptotic dynamical values, in particular the energy-density, 
coincide with the static ones though, the particle is not visiting any equilibrium state. 

In Fig. 1. we confront the analytical predictions (broken curves) to the numerical solution 
to the dynamical equations plotting m(t,t') vs. b(t,t') for t = 200,300,400. 7 = .5, T = 0.2, 
fi = 0.2 and 6 = 2.5 in the figure, q and the FDT relaxation converge rather quickly to the 
predicted value and straight line behaviour. However, the aging regime is still very far from 
the asymptotic analytical curve (for this choice of parameters the curve for the aging regime 
is unfortunately very similar to a straight line). The convergence for the long-range model 
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is slower than for the short-range model (see Fig. 2 and the discussion in Section VI. G). 
We observe that the evolution of the numerical curves as t increases is, though exceedingly 
slow, towards the dotted line. 

This family of long-range models behaves in every respect as the Sherrington-Kirkpatrick 
model of a spin-glass 11 . 



A. The phase diagram 



The discussion of the phase diagram for power law models with 7 < 1 is very different 
from the short range case. The function q 2 V"(q) is increasing monotonously thus there is 
no T^K The critical temperature T c (ji) is determined by: the condition b = q that implies 
V"{2T^ltj) = -// 2 /4 (see Eq.(VI.B.4). We can check that below this line \X[b]\ < 1/2T 



since: 



-2Tx(q) = - q - = < 1 (vii.i; 

W 2 V"(q) 2 9 + q y 



The dynamical phase diagram is thus identical to the statics. 
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VIII. SUMMARY AND CONCLUSIONS 



We have studied analytically the out of equilibrium dynamics of a particle moving in 
a general random potential. We have described in more detail three particular potential 
correlations associated to three well-known models: a power-law correlation, the Random 
Sine-Gordon model and the p-spin spherical spin-glass model. 

The model with power-type long-range correlation (model (1.1)) is solved by an ultra- 
metric ansatz for the triangular relations. In this case the asymptotic values of the one-time 
parameters such as the energy density, the equivalent of the Edwards- Anderson parameter 
q, the maximum displacement b 0} etc. coincide with the equilibrium values obtained with 
the full replica symmetry breaking scheme. Though the system never reaches equilibrium 
these quantities approach the equilibrium values asymptotically. This has been also shown 
to hold for the Sherrington-Kirkpatrick model and it is expected to hold in models without 
a threshold in the TAP energy landscape 10-12 . This seems to be a characteristic of models 
that are solved statically with a full replica symmetric scheme. Our results also agree with 
the numerical observations of Ref. [12]. 

Model (1.1) with short-range correlations is solved with a one blob ansatz for the tri- 
angular relations. This means that in the aging regime the function X[b] measuring the 
departure from the FDT theorem is simply a constant and that the dynamical equation de- 
termines the function / relating any three correlation functions. We have found that in most 
of the models we consider the inverse function / presents an irregularity at the begining of 
the aging regime. This irregularity is a new and interesting feature that, although generic, 
has not been previously observed in the off-equilibrium dynamics of mean-field models (e.g., 
it does not appear in the p-spin spherical model). A related irregularity, at the beginning 
of the out of equilibrium relaxations, has been predicted to occur in a model of traps with 
a broad distribution of trapping times by Bouchaud. This allows for a comparison of the 
mean-field glassy dynamics with the trap description of aging. 

54 



The one step trap model for the out of equilibrium relaxation of glassy systems 20 predicts 
that the integrated response of a system to a perturbation applied during the interval [0,1^] 
(the thermoremanent magnetisation for spin systems) has a very sharp decay to q in a 
microscopic characteristic time. Afterwards the aging regime sets on and the waiting time t w 
is its only characteristic time. The decay is then described by m(t, t w ) ~ (l — q)g(t)-\-qf(tjt w ) 
where the first term is the FDT relaxation and the second one is the aging relaxation. The 
aging part is characterised by two exponents 

f (f) ~ 1 ~ if) 1 '" 1 for *~*" f(f)~(rY' for *»*" ( VIIL1 ) 

Our results for short-range correlated potentials show that the mean-field dynamics of a 
particle in a random potential is also represented by a law of this type, at least in the 
beginning of the aging regime. Indeed when (t — t 1 ) jt 1 = 0(1) we find m(t, t') = x(b(t, t') — b ) 
which at the begining of this regime reads 



m(t, t)~x \ J + x ((l- b o)- 



(VIII.2) 



Thus, the exponent x\ is related to our exponent a. 

At the end of the aging regime, i.e. when t' jt ~ 0, we can compare with the behaviour 
obtained in the limit of a massless particle: 

&( M ') ~ 9 + In , (VIII.3) 

which naturally leads to logarithmic behaviour rather than power law at variance with 
(VIII.l). 



The RSGM model with non-zero mass belongs to the family of models solved by a one 
blob ansatz with irregularities. The RSGM in the limit of a vanishing mass as well as the 
p-spin spherical spin-glass are particularly simple; they can be studied in an independent 
way and they serve as checks to the method to obtain the exponent characterising the decay 
at the beginning of the aging regime. 
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As regards the case of a general correlation V we have obtained the conditions that V 
must satisfy in order to be solved by each particular ansatz. 

In many occasions it has been pointed out 40 ' 30 that the mean-field spin-glass dynamical 
equations, with the assumption of the relaxation being at equilibrium, are very similar to 
Gotze's mode coupling equations for a phenomenological description of the glass transition. 
The mode-coupling equations for glasses involve only the density-density correlation function 
and are time homogeneous by construction. Recently a very interesting first generalisation of 
the mode-coupling approach to account for non-equilibrium phenomena has been proposed 41 . 
The mean-field dynamical equations for the correlation and response function (without any 
time homogenous assumption) of the p-spin spherical spin-glass are indeed the dynamical 
equations arising from the mode-coupling treatment of a model without explicit quench 
disorder - the Amit-Roginsky <f> 3 model 41 . 

The equations we have studied in this paper constitute an enlargement of the class of 
possible extensions of the mode-coupling equations for glasses. Indeed, one can show that 
some of the mode-coupling equations of Gotze are those associated to the high-temperature 
phases of the models we consider here with particular choices of the random potential. 
Though in these models there is explicit quench disorder it is however interesting that the 
same equations appear. The further advantage of these equations is that they imply a more 
complicated behaviour for their solutions as compared to the solution of the p-spin spherical 
spin-glass, for instance, they have singularities at the begining of the aging regime that 
may be relevant experimentally for glasses 42 This issue will be discussed in more detail in a 
separate publication 43 . 

Finally this work opens the way to study other disordered models in finite dimension. 
We expect some features found here to extend to higher dimensions, such as aging effects 
with irregularity in the scaling function. Recent results on instabilities of renormalization 
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group flows in several models, such as the Sine-Gordon model 6 , suggest that aging effects 
will be present. 
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Appendix A 



In this Appendix we extend the Gaussian Variational Method (GVM) to the dynamics 
of disordered systems. This method was previously used to study the statics 5 . Here we 
extend it by applying a Gaussian decoupling to the exact equations of motion obtained from 
the Martin-Siggia-Rose generating functional. We derive in this way a set of dynamical 
equations associated to the model defined by the Hamiltonian (1.1) with a general Gaussian 
random potential V. 

We start by defining 

-S(x - x>) A(</> - </>') = (V(4>, x) V(4>', x')) . (A.l) 

The standard Martin-Siggia-Rose action entering the dynamical partition function = 
J T>cj) T>(f>exp( — S) associated to a Langevin process is 

r N 

S= dxdt J2[~T {ij> a {x,t)) 2 + ij) a {x,t) (<9 t -V 2 + /i) </> a (x,t)] 



r N r 1 * 

+ J dxdtdt 1 J2 "2 i^&t) Kp(<f>(x,t) ~ #M')) if 



af3 

s2 I £Ma.C±ll3 



(A.2) 



with A" /3 (</> — </)') = — (6 / 6 (f) a 6 (f) ) A((f) — </>'). The indices a, {3 label the components of the 
A^-dimensional vectors </>, icj). 

The dynamical equations for R a/3 (x, t; x' , f ) = (cf) a (x,t) i(f)P(x' ,t')), t > t', and 
C af3 (x,t;x',t') ee (<f> a (x,t) <p p (x' } t')) follow from 

/ i}P(x', t') S f ) = 6 aP 6 D (x - x') 6(t - t') , (A.3) 

\ oi<ft a (x,t) I 

(*"'■'• '''^H- (A - 4) 

respectively. These are exact equations of motion which can be derived by standard 
methods 44 . The brackets here and in what follows denote a mean over the fields weighted 
with the dynamical effective action. 
The first equation reads 
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S aP 6 D (x - x') 6(t - f) = (d t - V 2 + //) R afj (x,t; x' , f) 

-£ / dt"(z^(x',t') A'^OM) - if>(x,n) ) 

7 J 

(A.5) 

where we have used the fact that (icj)icj)) must be zero to preserve causality. The equation 
for C reads 

2T R af) {x\ t'; x, t) = {d t - V 2 + //) C af) {x, t; x', f) 

"E / dt"(r(x',t') A'^(x,t) - tW)) WW)) ) ■ (A.6) 

7 J 

The last term in the r.h.s. of the above equations can be computed using a Gaussian 
approximation, i.e. assuming that the fields (j) and icj) have a Gaussian distribution. Using 
the Ito prescription that states R(x } t; x } t) = 0, and the following rule for averaging any set 
of Gaussian variables j^} : (t^i* 1 (<£>)) = Y^ji^pi^Pj) the equation for R becomes 

6 af3 S D (x - x') S(t - t') = (d t - V 2 + //) R a ^(x,t; x' , t') 




x{A% c (</> -</>')). (A.7) 
Similarly, the equation for C reads: 

2T B aP (x',t';x,t)= (<9 t -V 2 + //) C aP (x, t; x\ f) - £ / dt"R^(x' ,t';x,t") (A'^ (<f> - <f>')) 

7 J 

~J2 J dt"R s ^(x, t; x, t") (C ta (x,t; x', f) - C ta (x,t"; x' , t')) 

x ( A< 4 > £ (</> -</>')). (A.8) 

where (..) denotes an average over a Gaussian distribution. These are the general dynamical 
equations obtained from the GVM. We now specialize to the case of O(N) symmetry, i.e. 
isotropy, which also includes N = 1 (RSGM) as a trivial case, and show that with the 
definition of V appropriate for studying the large N limit, the same function V has to be 
used in the statics and the dynamics. 
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Defining the two functions: 

-N V (^pj = A(</>) _JVV^0=(A(0)) (A.9) 

Using isotropy the equations (A. 7), (A. 8) contain only A"^ and A" ial3 which are related to 
V through: 

1 f d N q <1 2( P 2 
{A aP (<p)) = —6 aP J -^-jf A(q)q 2 exp(-^^) 

(E^(V))=4^V"(^) (A.10) 

7 

Replacing these expressions in (A. 7), (A. 8) one recovers the dynamical equations (2.1), (2. 2) 
with V substituted for V. In the limit of N — > oo these two functions become identical 5 . 

One can also apply the Gaussian decoupling in the exact equations of motion for the 
static replica theory. This provides an alternative way of obtaining the approximated saddle 
point equations of Ref. [5]. The exact saddle point equation is: 

n c I n \ 

(v 2 + ^)(Mx)Mo)) - (M*) E ( E V (M°) - Mo))j ) (A.n) 

Applying now the Gaussian decoupling one recovers the saddle point equations of Ref. [5]. 
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Appendix B 



In this Appendix we describe the derivation of the long-time dynamical equations, i.e. 
for times such that t > t' 3> 1 and it — t')/t' = 1. 
Neglecting the time-derivative, the r-Eq. reads: 

~ ~ r(t, f) (-n + 4 j*~ ds V"(B(t, s)) R(t, s) + 4 J* ds V"{b F {t - s))r F {t - s^j 
-4 T + ds V"(b(t, s)) r(t, s) r F (s - t') - 4 /*" ds V"(b(t, s)) r(t, s) r(s, t') 



t' Jt' + 

t 



4 / ds V"(b F (t - s j) r F (t - s) r(s, f) . (B.l) 



We have explicitly separated in the integrals the FDT regimes (integrals symbolically de- 
noted ds ) from the widely separated time regimes. In the last term of Eq. (B.l) s varies 
from t- to t and it is very far away from t'. The assumption is that in that case the functions 
vary very slowly (we have already neglected the time derivative in the l.h.s.). Hence r(s,t') 
is almost constant in this interval and it can be approximated by r(t, t'). The last term then 
cancels the third one. In addition, using the same approximation and the FDT relations 
(III. 5) one computes explicitly the integrals involving the FDT parts: 



-A ds V"(b(t, s)) r(t, s) r F (s - t') ~ -AV"(b(t, t'))r(t, t') / dr r F (r) 
Jt' Jo 

~ 2 V(HM')MM') 



T 

4 jf* ds V"(b f (t - s)) r F (t ~s) = ~l (V'(0) - V(q)) (B.2) 

using b F (oo) = q and b F (0) = 0. Finally, the second term inside the parenthesis involves 
both short and long times since the integral goes from the initial time, that is strictly zero 
to t- that is large and belongs to the asymptotic regime. One should then separate the 
contribution from finite times to that of long times 

f~ ds V"(B(t, s)) R(t, s) ~ /° + ds V"(b s (t, s)) r s (t, s) + f~ ds V"(b(t, s)) r(t, s) . (B.3) 
Jo Jo Jo + 
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We recall the weak long-term memory hypothesis and assume that the finite times do not 
contribute to the long-time dynamics and that the system forgets at large-times what hap- 
pened in the very short times after the initial time. 

Then the first term in the r.h.s. of Eq.(B.3) is neglected and one finally obtains the 
equation for the slow part of the response function: 

= r(t, f) + 4 J* ds V"(b(t, s))r(t, s) - ^V"(b(t, *'))) 

-4 /* ds V"(b(t, s)) r(t, s)r(s, t') (B.4) 
Jt' 

The same procedure can be carried through for b(t,t'). First one can use that at large 
t,t', C(t,t) ~ C(t',t') ~ cp(0) = q. Then, the integrals which appear in the 6-Eq. and 
contain FDT pieces are: 

2 J* ds V'(b F (t - s)) r F (t ~s) = ^- (V(0) - V(q)) 
-2 £' ds V'(b(t, 0) Mt' -s) = ^ V'(b(t, 0) 
2 J* ds V"(b F (t - s)) r F (t - s) b F (t - s) = | V'(q) + 1 (V(0) - V(q)) 
2 J* ds V"(b F (t - s)) r F (t - s) (b(t, f) - b(s, t')) ~ (B.5) 

One also needs to compute: 

-2 /* + ds V"(b(t, s)) r(t, s) b F (s, t') ~ -2V"(b(t, t')) r(t, t') f + ds b F (s, t') . (B.6) 
Jt'_ Jt'_ 

One can see that this integral is subdominant and gives a vanishing contribution to the 
b(t,t') equation in the large t limit. Since r(t,t') itself is of order 1/t (see below) it results 
that the FDT part of this integral is of order q(t' + — t'_)/t. The rapidly varying part of the 
total integral is thus vanishingly small compared to the total integral which is dominated 
by the slowly varying part. 

Again, this equation could have contributions from the initial times < s < + . We 
assume they are subleading, i.e. 
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~ / + ds V'(B(t, sj) (R(t, s) - R(t', sj) 
Jo 

~ / + ds V"{B(t, s)) R(t, s) (B(t, s) + B(t, f) - B(s, t')) . (B.7) 
Jo 

One finally obtains the equation for the slow part of the displacement correlation function: 

° = (" 2 + 2 fo dS V ' W ' S)) Kt ' S) ) b{t ' f ^ + T + ^ V '^ ~ V{h ^ *')) 

+2 /* ds V'(b(t, s)) r(t, s) - 2 /* ds V'(b(t, s)) r(t', s) 
Jo Jo 

+2 /* ds V"(b(t, s)) r(t, s)b(t, s) - 2 /* ds V"(b(t, s)) r(t, s)b(t', s) 
Jo Jo 

-2 /* ds V"(b(t, s)) r(t, s)b(s, f) . (B.8) 
Jt' 

From (2.3) one finds, using the same decomposition of the integrals and formulae (B.5) 
to evaluate the FDT parts, the large-time dynamical equation for C(t,t) and C(t,t'). 
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Appendix C 



Let us also indicate the derivation of another (equivalent) form for the equation for 
R(t,t') and r(t,t'). Starting from (2.1) one defines: 

F(t,t') = - f dsR(t,s) (C.l) 
Jt' 

i.e., such that R(t,t') = dF g t f ^ ■ One obtains: 

integrating over t' and using F(t',t') = one obtains: 

at Jo os 

-4 £ ds V"(B(t, s)) F(s, 0^7^ (C3) 

The integration constant has been fixed using the limit t — >■ t' and: 

v dF(t,t') f' dF{t,s) x 
hm )-^- = -bm (R (t,t + / ds 

f^t dt f^t K v ' ; J t > ds ' 

= -1 -lim(r F (0) -r F (e)) = -1 (C.4) 



Note that at this stage equation (C.3) is exact and provides another convenient form to 
equation (2.1). 
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Appendix D 



In this Appendix we obtain the recursion relations for the function B(u) suitable for 
a numerical solution. Defining v = u a and B(u) = vg(v) } one introduces g(v) n = 
Em=oi'( n ) m ) t ' m an d g( n 7®) = 1- Equation (VI. E. 6) then leads to: 

oo 

= cnv n g(v) n + Xc n . lV n H(n, v) (D.l) 

n=2 

where c n = V^ n+1 \q)/nl, X = —2Tx/q and 



oo 



1 d 



H(n,v) = h(n,m)v m = — — uv n / dxx {n -^ la {\ - x) a givx") 71 ' 1 g{v{\ - x) a ) . 

OU Gj%Xi iJ 

(D.2) 

After some algebra, we obtain two recursion relations for the coefficients g(n, m) and h{n, m): 

m 

g(n,m) = Y, 9(n-l,j)g(l,m-j) (D.3) 

3=0 

%,m = J2 9(n - l,j)g(l,m-j) 1 + /3 1 + J -, 1 + ^ , D.4 

„ a a a 

where f3(x, y) = T(x)T(y) /T(x + y). At the iteration n the coefficients g(i,j) and h(i,j) with 
i = 3, . . . , n, j = n — i = 0, . . . , n — 3 are first directly obtained from the recursion relations 
(D.3) and (D.4). The coefficients g(l,n - 2), g(2,n - 2) (and h(l,n - 2), h(2,n - 2)) are 
then obtained from the equation 

n 

2g(l, 0)g(l, n - 2) (c(2) + Xc(l)/3(1, n - 1)) = - ^ (c(%(&, n - k) + Xc(& - n - &)) 

fc=3 

ri — 3 

- X) flf(l, %(1, n - 2-k) (c(2) + Xc(l)/3(1 + jfc, n - 1 - &)) (D.5) 

it=i 

with 

n — 3 

5 (2, n - 2) = 20(1, 0)flf(l, ra - 2) + %(1, rc - 2 - &) (D.6) 

1 

n — 3 

/i(2,n - 2) = 2flf(l,0)flf(l,n - 2)/3(l,n - 1) + ^ #(1, %(1, n - 2 - k)/3(l + k,n - 1 - k) . (D.7) 

i 

where f3(x, y) = (1 + (x + y)/a)/3(l + x/a, 1 + j//a). 
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